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Abstract. 

We derive an explicit simple formula for expectations of all Schur functions in the 
real Ginibre ensemble. It is a positive integer for all entries of the partition even and 
zero otherwise. The result can be used to determine the average of any analytic series 
of elementary symmetric functions by Schur function expansion. 



PACS numbers: 02.10Yn, 02.50.-r, 05.40.-a, 75.10.Nr 
1. Introduction 

Real asymmetric random matrices, and in particular statistics of their (complex) 
eigenvalues, have many interesting applications in modelling of a wide range of physical 
phenomena. They appeared in the studies of stability of complex biological networks 
PQ, dynamics of neural networks with asymmetric synaptic couplings [2], directed 
quantum chaos in randomly pinned superconducting vortices [3], delayed time series 
in financial markets [I], quantum chromo dynamics with a real representation of the 
Dirac operator [5], or random quantum operations in quantum information theory [6]. 

All statistical information about eigenvalues of random matrices is encoded in 
their joint probability density function (pdf). For real asymmetric random matrices 
the joint pdf of eigenvalues was derived in [7] almost 20 years ago. It gave insights 
about the generic properties of the eigenvalue distribution, like the level repulsion, but 
had a complex structure due to singularities which made the precise statistical inference 
about eigenvalues difficult at that time. It suffices to say that the eigenvalue counting 
measures for complex [8] and real eigenvalues [9] were initially obtained without the 
aid of knowing the joint pdf. Very recently the interest in the spectral properties of 
real asymmetric random matrices has been revived [10] and the Pfaffian representation 
of n-point eigenvalue correlation functions has been derived directly from the joint 
pdf by various methods [HJ [121 021 EE]- This progress made it possible to compute 
the bulk and edge scaling limits of the eigenvalue correlation functions [15] and, in 
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particular, to recover the results of [SJ [§] about the eigenvalue counting measures in the 
limit of infinite matrix dimension [T2| [T3] . Interestingly, the kernel that determines the 
eigenvalue correlation functions can be written in terms of averages of the characteristic 
polynomial of random matrices of lower dimension, much in the spirit of [H[9], see also 
[T6] for complex matrices, and hence can be computed in a very simple way |17J . One 
intriguing open question is the relation of real random matrices to integrability theory 
such as r-functions which are generating functions for the cumulants of power sums of 
eigenvalues. 

The eigenvalues of the real asymmetric matrices are either real or pairwise 
complex conjugate. This implies a complicated and singular structure of the eigenvalue 
correlation functions. It seems that in many aspects the elementary symmetric functions 
of eigenvalues, e n , or their first N power sums, t n , form a more convenient set of 
independent real variables. The Jacobian of the transformation from this set to 
the eigenvalues is just the Vandermonde determinant appearing in the joint pdf of 
eigenvalues [?]. Since symmetric functions of eigenvalues can be expanded in products 
of e n 's (or t n 's), averages of symmetric functions of eigenvalues can be written in terms 
of the joint moments of elementary symmetric functions (or power sums). Equivalent 
to all these moments is the set of all averages of Schur functions which form another 
basis in the space of symmetric functions. As Schur functions are determinants, divided 
by the Vandermonde their average can again be reduced to a Pfaffian. This Pffafian 
can be evaluated explicitly as we will show in this paper. The result is very simple. If 
the partition characterizing the Schur function is odd then the Schur function average 
is zero. For even partitions the Schur function average is a positive integer expressible 
in terms of Gamma functions. This gives an alternative method of evaluating averages 
of symmetric functions of eigenvalues by the way of Schur function expansions. 

We will treat both even and odd matrix dimensions iV in a coherent way. We will 
reveal a further interesting structure of the theory. The basic building block is a skew 
symmetric kernel K,n(z\, z 2 ), which is a polynomial of degree N — 1 in both eigenvalues 
zi, Z2, thus building a skew symmetric from with a matrix A" 1 . If one knows this matrix 
A" 1 in all even dimensions, one is able to determine all correlation functions and all 
moments. Since A -1 has the same form in all dimensions, which just can be considered 
as a submatrix of one infinite dimensional matrix, the averages of general Schur functions 
are given by Pfaffians of finite submatrices of this infinite dimensional matrix. For odd 
dimension one has to modify this procedure slightly. 

2. The real Ginibre ensemble 

In this paper we consider real random N x N matrices H with probability measure 




(1) 
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This ensemble of random matrices is known under the name of real Ginibre ensemble 
[T8] . The matrix distribution (JTJ induces a probability distribution on the eigenvalues 
of H which is known to be of the form [3, M, EE] 

dfi(zi, z 2 ,..., z N ) = C N ■ Y[(zi - Zj) ■ Yl f( z k) ■dz 1 ... dz N , (2) 

i<j k 

with positive function f{z) = f{z) > 0, f 2 {z) = e~ Re(z2 ) erfc(| lmz\y/2). The eigenvalues 
Zj of H are real or pairwise complex conjugate and are ordered in such a way that dfi > 0. 
Introducing the skew symmetric weight function (z = x + iy, d 2 z = dxdy) 

F(zi, z 2 ) = f(z 1 )f(z 2 )(2i5 2 (z 1 - z 2 )sgn(y 1 - y 2 ) + 5(yi)5(?/2)sgn(a;2 - x{)) 

with 5 2 {zi — z 2 ) = 8{x\ — x 2 )6(yi + y 2 ), all correlation functions are determined with 
the help of the skew symmetric kernel 

N 

JCn(zi,z 2 ) = ^zt'zif 1 (3) 

k,l=l 

with the skew symmetric matrix 

A kl = J d 2 Zl J d 2 z 2 T{z^z 2 )z\- l z l 2 - x . (4) 

The eigenvalue correlations are Pfaffians with combinations of /Cat and JF as entries 
[19]. Comparing the 1-point density R\(z) = J d 2 uJ-'{z 1 u)lC^{u 1 z) with Edelman's 
expression [8] for the density of complex eigenvalues, 

one finds 

vi \ z i~ z 2 sr^ \z\z 2 ) , 

v n=0 

This result can independently be found by a simple supersymmetric (Grassmannian) 
calculation of the average of a product of two characteristic polynomials. This 
calculation can be extended to partly symmetric or chiral counterparts [17J. 
With the help of the duplication formula for the Gamma function 

2v / 27r(AT - 1) = 2( iV - 1 )/ 2 r((iV - l)/2) • 2 N ' 2 Y{N/2) (6) 

([3]) and (jSJ) imply that 

A u = o-k^kiaY 1 (7) 

with cifc = 2 k / 2 T(k/2) and e^i being the tridiagonal skew symmetric matrix with —1 in 
the upper diagonal. The structure of the Vandermonde determinant in (j2j) and (or) the 
results of paper [19] imply that this expression can be used for odd dimension too, if we 
extend the matrix A^ 1 to dimension N + 1 by putting otv+i = 1 in (|7j). In this case the 
skew symmetric matrix e can be inverted and the inverse yields (Hj) from (JTj). However, 
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in that case formula ([3]) holds with the sum over k, I running only from 1 to N. The 
normalization constant C/v in (T5]) is given by 

\/C N = Pfaff(A) (8) 

that is the Pfaffian of the skew symmetric matrix A. 



3. Schur function averages 

The description of the system by the complex eigenvalues might not be the most 
convenient one. A more suitable set of independent real variables is the set of power 
sums 

TV 

t n = Tr(H) n = E Z i ( 9 ) 

i=l 

or the set of elementary symmetric functions ei, . . . , ejv which are the coefficients in the 
expansion of the characteristic polynomial in powers of x, 

N 

det{l + xH) = ^x n e n {z 1 ,...,z N ) . (10) 

n=0 

Sometimes complete symmetric functions h n come in handy These are the coefficients 
in the expansion of the reciprocal characteristic polynomial in powers of x, 

oo 

1/ det(l -xH) = J2 ^M^i, ...,z N ). (11) 

n=0 

Homogeneous symmetric polynomials in Zi of degree > 1 can be expressed as sums of 
products of power sums t n , or equally as sums of products of e n 's or, equally, of h n 's. 
Hence, when calculating averages of the symmetric polynomials in eigenvalues of H, 
or analytic functions, one is led to the problem of finding the joint moments of the 
eigenvalues. A full description of all symmetric moments of eigenvalues can be given 
by the so-called Schur polynomials or Schur functions. These functions are indexed 
by sequences, called partitions, A = (Ai,...,A/v) of ordered non- negative integers 
Ai > A 2 > • • • > Ajv > 0, called parts, and can be conveniently written as the ratio of 
two determinants 

a x (z u ...,z N ) =det(^-" +A ")/det(^- n ) . (12) 

We will follow the convention of not showing zero parts, i.e. 

(Ai, A2, • • • , A/, 0, . . . , 0) = (Ai, A2, • • • , A/) 

and when referring to specific partitions will adopt the notation A = (k ak , . . . , 2 a2 , l ai ) 
to indicate that there are precisely a\ of l's, a 2 of 2's, etc., among the parts of partition 
A. 

It can be verified directly from ( |T2l) that if partition A consists of one part, 
A = (n), then the corresponding Schur function is just the complete symmetric function, 
C(n) = hn and if partition A has only zeros or ones among its parts, A = (l n ) then the 



Schur function averages for the real Ginibre ensemble 



5 



corresponding Schur function is just the elementary symmetric function, o~(rn\ = e n . 
Also, by convention, a\ = is the number of non-zero parts of A is greater than the 
number of indeterminants Zj. 

If A is a partition of n, i.e. A x + . . . + Ajv = n, then the Schur function o~\ is 
a homogeneous symmetric polynomial of degree n. In fact the set of all such Schur 
functions form a basis in the space of homogeneous symmetric polynomials of degree n. 
In particular, for n > 1, see, e.g., [20J, 

n 

t n = , (13) 

fc=l 

where the sum is over all hook partitions A = (k, l n_fc ) of n. Vice versa, as already 
mentioned above, every Schur function, as any other symmetric polynomial, can be 
expressed in terms of products of power sums. One consequence of this is that the 
Schur function of matrix argument defined by o~\(H) = a\(zi, . . . , zn) with Zj being the 
eigenvalues of H, is a polynomial in the matrix entries of H and o~\(H) = o-\(THT~ r ) 
for any non-degenerate matrix T. 

The Schur functions are the characters of irreducible representations of the unitary 
group and as such are orthogonal with respect to integration over the Haar measure. 
This property is quite useful for finding the coefficients in Schur function expansions. 
For example, expanding the product of characteristic polynomials in Schur functions, 

n N 

\ | + XjZ k ) = ^ c \( x l, • • • ; Xn)<7x(zi, ■••,Zn), 

j=l k=l A 

one can find the "Fourier coefficients" c\(xi, . . . ,x n ) by integration over the unitary 
group 

/n 
Y[det(I + Xj U) ax(U)dU, (14) 
i=i 

where the integration is over the unitary group U(N) with respect to the normalized 
Haar measure dll. The group integral in (fl4"l) can be easily evaluated yielding 
ca = det(e\ k -f-+j) which, by the Jacobi-Trudi formula is again a Schur function, 
det(eA fc _fc+j) = o"a', with A' being the partition conjugate to Thus one arrives 
at the Schur function expansion 

n N 

a y (x 1 ,...,x n )a x {H). (15) 

j=l k=l A 

This expansion is well known under the name of dual Cauchy identity and can be derived 
by purely algebraic means, see e.g. [21]. Other known Schur function expansions can 
be obtained via ( TT4l) . see, e.g., Appendix in [22] . 

The Schur function averages over the real Ginibre ensemble ([1]) can be obtained 
by the same reasoning as in (J^J) - §E§ or following the calculations in [19|. The only 

| The conjugate to partition A = (Ai, A2, . . .) is the partition A' = (X[, X' 2 , . . .) with A^- being the number 
of Afe that are greater than or equal to j. 
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difference is that the Vandermonde det(z^ n ) has to be replaced by det(^ n+A «). The 
result is again a Pfaffian, 

(a x (H)) N oc Pfaff J d 2 ' Zl d 2 'z 2 T \z u z 2 ) z$ r ~ n+Xn z* r ~ m+Xm (16) 

where the brackets mean average over the real Ginibre ensemble (CQ) and n,m = 
1,2... AT. The normalization constant can be restored requiring that (o~\(H))n = 1 
for the empty partition A = (0). By making use of (TjJ and (JTj), one can also write 

(a x (H)) N oc Pfaff A A r_ n+An +i,Ar_ m+ A m +i • (17) 

The indices N — n + X n + 1 form an increasing sequence in the opposite direction: 

1 + \ N < 2 + Ajv-i < ... < N + Ai . (18) 

Here it does not matter that N — m + X m + 1 can be > A" because it turns out that 
A^ m has the same form in all dimensions. The reason is that {e)~ l m can be obtained 
from any M dimensional matrix e n ^ m with even M > n,m. This is due to the special 
tridiagonal structure, as one can easily check (see also the following). Therefore, and 
this is the simple idea, we can immediately calculate all 



.4 



l:m 



a l ( e )hm a m 



(19) 

which then is valid also for any submatrix centered along the diagonal. Again for odd 
dimension A^ we obtain e _1 from a higher even dimension say M cutting all rows and 
columns with number not equal N — n + X n + lorM and put au = 1- This is meant-and 
correspondingly in the even A^ case- if we write e~^_ n+Xn+1 N-m+\ m - Thus 

(ax{H)) N oc Pfaff(a Ar _ n+An+ ie^ 1 „ n+An+1 Ar _ m+Am+1 a A r_ m+Am+ i) . (20) 
To calculate the Pfaffian, let us define it by a Grassmann integral (or Berezin integral) 

Pfaff(v4 ijm ) = J dx\... dxN exp(-^ XiA,mXm) (21) 



Lm 



with anticommuting variables Xi- Then the Berezinian (Jacobian) for the transformation 
Xi ~^ Xi/ a i is Yli a i = Ylb=i a i ( n °te that it appears with the opposite exponent as 
compared to commuting variables). Thus we obtain 

TV 

(cx{H))n OC a A r -n+A„+l)Pfaff ( e N-n+\ n +l,N-m+\m+l) ' 

(22) 

n=l 

For calculating the matrix (ejjTj) we may choose the smallest even dimension M, that 
contains all N - n + A n + 1 : M >N + X 1 + 1. 

Let us illustrate the matrix e, which we always use in even dimension 



/ o 


-1 
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-1 
















+1 





-1 
















+1 
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■•-/ 



(23) 
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It is easy to see that Pfaff(— e) = J dxi---dxNi 
corresponding inverse e _1 has the form 



-1 


-1 



-1 

V : 



+1 














+1 



+1 















+1 



+1 



+1 





X1X2X-X3X4) 

\ 



1 The 



(24) 



/ 



Again we see that Pfaff(e _1 ) = f dxi ■ ■ ■ oIxn{.— XiX2){~ X3X4) • • • — 1- Now let us 
calculate the Pfaffian in (1221) 
simpler (remember (PfafL4) 2 



We first try to calculate the determinant which is 
det A) , then we determine the sign. To calculate 
subdeterminants of e _1 we consider the generating function 

D N = det(x fc 5 M + (e -1 )*,,) . (25) 

From this we can calculate subdeterminants by differentiating with respect to x k 
to cut row and column number k and then put x = 0. Now we transform using 
det e _1 = det e = 1 



N 1 
D N = (TT x k ) det(— S k 1 + e k>i 



(26) 



k=l 



which now we may generalize to odd N and deduce due to the simple form of e 

D N = Djv-i + x N x N -iD N _ 2 (27) 

with Di — 1, D 2 = 1 + x\Xi- As a result 

D N = 1 + xix 2 + x 2 x s + . . . + £1X2X3X4 + . . . (28) 

is a sum of products of consecutive pairings with all coefficients equal to +1. This implies 
that we can cut only consecutive pairs. As a result: |Pfaff(e^ v 1 i n+A)i+1 N _ m+X = 1, 
if all X n even and = if one or more \ n odd. To determine the sign we observe 
that a permutation of rows and columns of the Pfaffian changes it by the sign of the 
permutation. Let us move in this way all rows and columns, which we don't want to cut, 
into the upper left corner of the matrix. This yields a sign (— 1)' A ' of the permutation 
(which is = +1 if |A| = Ai + A2 + • • • + Aat is even). Then we observe, that if we move 
away only consecutive pairs of rows and columns, the remaining submatrix of e" 1 in the 
upper left corner does not change its form. Thus its Pfaffian is +1. As the final result 
we obtain 



2 2 







N 

n 

n=l 



r((iV-7i + A w + l)/2) 
r((JV-n + l)/2) 



if all A„ are even 



otherwise. 



(29) 



Here |A| = Ai + . . . + Ajv and the upper limit N of the product can also be replaced 
by the length of the partition, which is the number of its nonzero parts A n . It has to 
be mentioned that the Schur function average (1291) . can be also obtained from zonal 
polynomials [21], however this derivation is rather indirect. 
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4. Conclusions 

We have calculated the average of Schur functions over the real Ginibre ensemble by 
using the Pfafhan representation fl20]) for the average and then evaluating the Pfaffian 
explicitly. It vanishes for odd partitions and is given by a simple expression in terms 
of Gamma functions ( 1291) for even partitions. This result might be useful for evaluating 
averages of symmetric polynomials and analytic functions in eigenvalues of real Gaussian 
matrices. For example, consider the power sums t n , n > 1. Their averages can easily be 
calculated by the way of Schur function expansion ffT3|) and our main result fl29|) . For 
odd n there are no hook partitions of n with all parts even and therefore the average of 
t n vanishes. Of course, this can be verified directly from the symmetries in the matrix 
distribution ([1]). If n is even then there is only one hook partition of n with all parts 
even. This partition is A = (n) . Hence, all terms but one on the right-hand side in (fl3l) 
vanish and one arrives at the relation 

(TrH n ) N = (a (n) (H)) N = (h n (H)) N 

where h n is the complete symmetric function of degree n. Hence, by ( 1291) . 



(Tr H^)n = 2 m r( r^ } m) = IK N + 2 ( m ~ j)l 
Similarly, on averaging (115j) with the help of (|29l) one obtains 



3=1 n/2>\ 1 >X 2 >...>X N >0 



where A = (A 1; ...,X N ),2X = (2A 1; . . . , 2X N ) and c A = 2^ U k T( r«N-k+W) 2) • If n = 2 
then the constraints n/2>Ai>A2>...>Ajv>0 imply that Xj = 1 or for each 
j. Correspondingly, A = (l fc ), 2A = (2 k ) and (2A)' = (k,k), where k is the number 
of non-zero Xj. Observing that o~^,k)( x i: x 2) — {xi%2) k an d cnk\ = N\/(N — k)\ one 
concludes that 

N ( \ k 

(det(l + Xl H) det(l + x 2 H)) N = N\ ^ ( % k)V (30) 

k=0 \ >' 

in agreement with calculations in [8] and [17]. This expansion can also be obtained via 
(flQj) directly from the symmetries in the ensemble distribution ([1]). One only needs to 
recall the fact that the elementary symmetric polynomial of degree n in the eigenvalues 
of H is the sum of all principal minors of H of order n. Then fl30l) follows from the 
invariance of the distribution ([T]) with respect to changing the sign of Hi m and, also, 
permutation of columns of H. 
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